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The  search  for  clean  and  renewable  sources  of  energy  represents  one  of  the  most  vital  challenges  facing 
us  today.  Solid  oxide  fuel  cells  (SOFCs)  are  among  the  most  promising  technologies  for  a  clean  and  secure 
energy  future  due  to  their  high  energy  efficiency  and  excellent  fuel  flexibility  (e.g.,  direct  utilization  of 
hydrocarbons  or  renewable  fuels).  To  make  SOFCs  economically  competitive,  however,  development  of 
new  materials  for  low-temperature  operation  is  essential.  Here  we  report  our  results  on  a  computational 
study  to  achieve  rational  design  of  SOFC  cathodes  with  fast  oxygen  reduction  kinetics  and  rapid  ionic 
transport.  Results  suggest  that  surface  catalytic  properties  are  strongly  correlated  with  the  bulk  transport 
properties  in  several  material  systems  with  the  formula  of  Lao.5Sro.5BO2.75  (where  B  =  Cr,  Mn,  Fe,  or  Co). 
The  predictions  seem  to  agree  qualitatively  with  available  experimental  results  on  these  materials.  This 
computational  screening  technique  may  guide  us  to  search  for  high-efficiency  cathode  materials  for  a 
new  generation  of  SOFCs. 

©  2009  Published  by  Elsevier  B.V. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFOs)  have  attracted  considerable  atten¬ 
tion  due  to  their  high  energy  efficiency  and  excellent  fuel  flexibility 
(e.g.,  the  potential  for  direct  utilization  of  hydrocarbons,  coal  gas, 
biomass,  and  other  renewable  fuels  [1,2]).  To  realize  these  unique 
advantages  over  other  types  of  fuel  cells,  however,  several  critical 
challenges  must  be  overcome,  including  anode  deactivation  by  cok¬ 
ing  and  contaminates  (e.g.,  sulfur)  poisoning,  cathode  degradation 
by  chromia  poisoning,  and  high  cost  for  broad  commercializa¬ 
tion  [3].  One  effective  approach  to  cost  reduction  is  to  reduce 
the  operating  temperature  of  SOFCs.  At  sufficiently  low  operat¬ 
ing  temperatures,  the  interconnect,  heat  exchangers,  and  other 
structural  components  can  be  fabricated  from  much  less  expen¬ 
sive  materials.  Further,  as  the  operating  temperature  is  reduced, 
many  technical  difficulties  will  find  easy  solutions,  system  reli¬ 
ability  and  operational  life  increase,  and  so  is  the  possibility  of 
using  SOFCs  for  a  wider  range  of  applications.  It  is  noted,  how¬ 
ever,  that  the  interfacial  polarization  resistances  increase  rapidly 
as  the  operating  temperature  is  reduced,  especially  the  resistances 
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to  oxygen  reduction  on  the  cathode  [4],  implying  that  the  cre¬ 
ation  of  novel  cathode  materials  with  high  catalytic  activity  toward 
oxygen  reduction  at  low  temperatures  represents  a  critical  step 
toward  the  development  of  low-cost  SOFCs  [5].  To  reduce  cathodic 
polarization,  mixed  ionic-electronic  conductors  (MIECs)  such  as 
(Lai_xSrx)(Bi_yB^)03_5  (B  or  B'  =  Mn,  Co,  or  Fe)  have  been  widely 
used  as  cathode  for  SOFCs.  It  has  been  reported  that  LaCo03 -based 
cathode  materials  exhibit  better  performance  than  LaMn03 -based 
ones  [6],  although  the  detailed  mechanisms  of  oxygen  reduction 
on  these  materials  are  still  not  clear.  Recently,  density  functional 
theory  (DFT)  calculations  have  been  used  to  examine  the  ener¬ 
getics  for  the  elementary  steps  involved  in  oxygen  reduction  on 
LaMn03 -based  cathodes  (LaMn03  and  Lao.5Sro.5MnO2.75)  [7-9], 
including  adsorption,  dissociation,  incorporation,  and  diffusion  to 
triple-phase  boundaries  (TPBs)  where  a  cathode,  an  electrolyte, 
and  oxygen  species  meet.  In  addition  to  the  kinetics  of  the  oxy¬ 
gen  reduction  reaction  on  cathode  surfaces,  transport  of  oxygen 
ions  to/away  from  the  surface  active  sites  plays  an  important  role 
in  rational  design  of  novel  cathode  materials  for  SOFCs  [  1 0,1 1  ].  Sec¬ 
ondary  ion  mass  spectrometry  (SIMS)  [10,12,13]  is  considered  one 
of  the  most  reliable  experimental  techniques  for  determination  of 
oxygen  ion  diffusivity  and  catalytic  properties  (e.g.,  surface  oxygen 
exchange  coefficient)  of  MIEC  cathode  materials.  Further,  a  series  of 
atomistic  simulations  [14-17]  have  also  been  used  to  predict  ionic 
transport  through  AB03-type  materials.  In  this  study,  we  report  our 
results  on  a  theoretical  prediction,  on  the  basis  of  first-principles 
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methods,  of  both  surface  reactions  and  bulk  diffusion  in  Sr-doped 
LaB03-based  cathode  materials  (B  =  Cr,  Mn,  Fe,  or  Co).  In  particular, 
we  have  investigated  how  ionic  diffusion  through  bulk  cathodes 
and  surface  interactions  between  oxygen  species  at  their  surfaces 
influence  the  efficiency  of  the  cathode  performance  in  an  SOFC,  pro¬ 
viding  important  insight  into  rational  design  of  high-performance 
cathode  materials  for  SOFCs. 

2.  Computational 

For  electronic-structure  calculations,  the  Vienna  ab  initio  simu¬ 
lation  package  (VASP)  [18,19]  was  used  with  projector  augmented 
wave  (PAW)  potentials  [20].  La,  Sr,  Cr,  Mn,  Fe,  Co,  and  O  atoms 
were  described  by  11  (5s25p65d16s2),  10  (4s24p65s2),  6  (4s13d5),  7 
(4s13d6),  8  (4s13d7),  9  (4s13d8),  and  6  (2s22p4)  valence  electrons, 
respectively,  while  Ce,  Pr,  Nd,  Pm,  and  Sm  atoms  were  described  by 
12  (5s25p64f15d16s2),  11  (5s25p65d16s2),  11  (5s25p65d16s2),  11 
(5s25p65d1 6s2 ),  1 1  (5s25p65d1 6s2),  and  1 0  (5s25p65d1 6s1 ),  respec¬ 
tively.  The  Perdew-Wang  (PW91 )  exchange-correlation  functional 
[21]  was  applied.  All  calculations  were  performed  using  the  elec¬ 
tron  spin  polarization  with  a  400  eV  kinetic  energy  cutoff  for  a  plane 
wave  basis  set  and  a  Monkhorst-Pack  [22]  mesh  with  (4x4x4)  k- 
points.  The  Jahn-Teller  distortion  may  not  be  critical  under  SOFC 
operating  conditions  (above  500  °C  at  ambient  air)  [23,24]  since 
LaMn03 -based  materials  have  a  cubic  crystal  structure.  In  addition, 
Predith  and  Ceder  [25]  successfully  demonstrated  oxygen  ionic 
transport  phenomena  through  cubic  perovskite-structure  AB03 
(A  =  Ba  or  Sr  and  B  =  Mn,  Fe,  Co  or  Ru)  using  DFT  calculations.  Thus, 
we  applied  the  Pm3m  crystal  structure  to  construct  La-based  sur¬ 
face  models  similar  to  previous  studies  [26-30]  on  LaMn03-based 
materials.  Regarding  the  magnetic  properties,  only  FM  states  were 


applied  in  this  study  because  ferromagnetic  (FM)  configurations 
are  more  stable  than  antiferromagnetic  (AFM)  configurations  with 
a  difference  in  energy  of  0.02-0.25  eV,  depending  on  B  cations. 

3.  Results  and  discussion 

3.1.  Prediction  of  bulk- diffusion  barriers 

The  validity  of  the  surface  models  used  in  this  study  was  con¬ 
firmed  by  comparing  the  predicted  diffusion  barriers  with  those 
determined  from  experimental  SIMS  results.  All  bulk-diffusion  cal¬ 
culations  were  carried  out  on  Sr-doped  LaB03(l  1  0)  since  it  is 
reported  [15]  that  ionic  diffusion  may  occur  most  likely  in  the  (110) 
direction.  Because  the  defect  structures  of  (Lai_xSrx)(Bi_yB^)03_5 
(B  or  B'  =  Mn,  Co,  or  Fe)  depend  strongly  on  temperature  and 
partial  pressure  of  oxygen  [31,32],  we  constructed  a  simplified 
surface  model  with  a  stoichiometry  of  Lao.5Sro.5BO2.75.  Four  sur¬ 
face  models,  Lao.5Sro.5CrO2.75  (LSCr50),  Lao.5Sro.5MnO2.75  (LSM50), 
Lao.5Sro.5FeO2.75  (LSF50),  and  Lao.5Sro.5CoO2.75  (LSC50),  were  con¬ 
structed  according  to  the  following  defect  reaction: 

2SrO  +  2La£a  +  Oq  ->  2Sr^a  +  V0##  +  La203, 

where  La£a,  Oq,  Sr^a,  and  V0*#  are  a  La  cation  at  a  regular  La  site,  an 
oxygen  ion  at  a  regular  oxygen  site  in  the  bulk  phase,  a  Sr  cation 
at  a  La  site,  and  an  oxygen  vacancy,  respectively.  Recently,  Porn- 
prasertshk  et  al.  [33]  reported  that  diffusion  barriers  of  oxygen 
ions  through  bulk  phases  are  influenced  by  locations  of  oxygen 
vacancies  and  neighboring  ions,  leading  to  enormous  number  of 
possible  configurations.  To  minimize  the  number  of  possible  sur¬ 
face  models  and  to  examine  the  effect  of  oxygen  vacancies  on  the 


Fig.  1.  Schematic  of  a  surface  model  for  the  examination  of  ionic  diffusion  through  bulk  phases,  (a)  The  Lao.sSro.sMnC^jsCl  1  0)  surface  model,  (b)  Top  and  side  views  of  a  bulk 
Lao.5Sr0.5Mn02.75(l  1  0)  structure  before  and  after  ionic  conduction.  The  supercell  used  in  this  study  contains  a  total  of  19  ions  (2  La,  2  Sr,  4  B,  and  110  ions)  and  an  oxygen 
vacancy.  V  and  0N  denote  an  oxygen  vacancy  and  the  nearest  neighboring  oxygen  to  V,  respectively,  (c)  Trajectory  of  oxygen  ion  conduction  through  Lao.5Sro.5Mn02.75(l  1  0). 
01  and  02  denote  initial  and  final  states  of  oxygen  ion  conduction.  For  clarity,  oxygen  ions  involving  bulk  diffusion  are  presented  in  yellow.  (For  interpretation  of  the 
references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 
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Table  1 

Predicted  properties  and  diffusion  barriers  of  oxygen  ionic  conduction  through  Lao.5Sro.5BO2.75  (B  =  Cr,  Mn,  Fe,  or  Co)  at  the  GGA-PW91  level. 


Lao.5Sro.5CrO2.75 

Lao.5  Sro.5  Mn02.75 

Lao.5Sro.5FeO2.75 

Lao.5  Sro.5  C0O2.75 

Bulk  diffusion 

Em 

0.75 

0.76 

0.42 

0.38 

E  O-vac 

3.56 

3.25 

2.41 

1.49 

Ea.calc. 

4.31 

4.01 

2.83 

1.87 

Fa,expt.a 

- 

3.76  ±0.33 

- 

1.41  ±0.02 

Surface  interactions 

AEads,  superoxo- 

-1.86 

-1.37 

-1.07 

-0.82 

A£ads,  peroxo- 

-2.03 

-1.20 

-0.78 

-0.09 

A£ads,  mean 

-1.95 

-1.28 

-0.92 

-0.45 

a  From  Refs.  [12,13]. 


02 -cathode  interactions,  we  located  an  oxygen  vacancy  on  the  top 
layer,  as  shown  in  Fig.  1  a.  Furthermore,  similar  to  the  previous  stud¬ 
ies  [33,34],  we  assumed  that  the  diffusion  barrier  heights  depend 
strongly  on  the  nearest  oxygen  ion  site  to  an  oxygen  vacancy.  The 
ionic  conductivity  (D)  can  be  predicted  using  transition  state  theory 
(TST)  [35]  according  to 

D=D°exp(-|?)’ 


where  D0,  Ea,  kB,  and  T denote  a  pre-exponential  factor,  an  diffusion 
barrier  in  eV  at  0  K,  Boltzmann’s  constant,  and  temperature  in  K, 
respectively.  Diffusion  barrier  Ea  can  further  be  described  [36]  by 

Ea  —  Em  +  Eo-vac  > 

where  Em  and  Eo-vac  are  a  migration  barrier  of  an  oxygen  ion 
through  bulk  phases  and  an  oxygen-vacancy  formation  energy, 
respectively.  On  the  basis  of  the  previous  work  for  predicting  the 
migration  barriers  for  perovskite-type  materials  by  means  of  atom- 


Fig.  2.  Illustration  of  surface  oxygen  species  on  the  LSM50(1  0  0)  surface.  The  colors  illustrating  atoms  are  the  same  as  those  used  in  Fig.  1.  (a)  A  representative  surface  model 
with  four  layers  for  the  interactions  between  oxygen  species  and  Lao.5Sro.5MnO2.75.  (b)  The  species  in  dashed  circles  are  surface  oxygen  species  formed  via  adsorption  or 
incorporation  into  oxygen  vacancies.  V  denotes  an  oxygen  vacancy,  (c)  Adsorption  energies  of  surface  oxygen  species  versus  B  cations.  Averaged  02  adsorption  corresponds  to 
averaged  adsorption  energies  of  superoxo-  and  peroxo-like  species.  “Dissociated”  denotes  the  adsorption  energies  of  two  dissociated  oxygen  species,  while  and  “Incorporated” 
is  that  one  of  dissociated  oxygen  species  is  incorporated  into  an  oxygen  vacancy,  (d)  A  correlation  between  surface  diffusion  barriers  of  adsorbed  0  species  on  B  cations 
( AEdjff)  and  averaged  O2  adsorption  energies  on  four  Sr-doped  LaB02.75  surfaces  ( AEadSj0 ). 


1444 


Y.  Choi  et  al.  /  Journal  of  Power  Sources  195  (2010)  1441-1445 


istic  modeling  techniques  [15-17],  we  predicted  barrier  heights  of 
oxygen  ion  migration  using  the  nudged  elastic  band  (NEB)  method 
[37].  Shown  in  Fig.  lb  are  the  initial  and  final  sites  of  oxygen  ions 
in  the  LSM50  surface  model,  clearly  illustrating  a  curved  trajec¬ 
tory  of  the  oxygen  vacancies  from  position  01  to  position  02.  The 
non-straight  trajectory  is  in  agreement  with  previous  studies  [14], 
where  they  predicted  a  0.86  eV  migration  barrier  of  oxygen  ions 
through  the  LaMn03  bulk  using  atomistic  modeling  techniques. 
However,  it  is  not  straightforward  to  compare  with  our  results 
because  dopant  Sr2+  cations  influence  the  barrier  height.  In  addi¬ 
tion,  as  Islam  discussed  [15],  because  the  predicted  value  from 
atomistic  modeling  is  from  only  oxygen-vacancy  migrations,  they 
are  much  lower  than  those  determined  from  SIMS  data— De  Souza 
and  Kilner  [12,13]  reported  that  an  experimental  activation  energy 
of  LSM50  from  the  SIMS  technique  is  3.76  ±  0.33  eV.  In  order  to 
predict  its  diffusion  barrier,  we  estimated  E0_Vac  (Lao.5Sro.5BO3  -► 
Lao.5Sro.5BO2. 75  +  V0-  +  (l/2)02(g))  according  to 

Lo-vac  —  LILao.5Sro.5BO2. 75]  +  ^E[ 02]  -EILao.5Sro.5BO3], 

where  Lao.5Sro.5BO3  and  Lao.5Sro.53MO2.75  denote  Sr-doped  La- 
based  cathodes  without  and  with  an  oxygen  vacancy  (V0— ), 
respectively  [38],  while  E[02]  is  the  predicted  electronic  energy  of 
triplet  02  in  a  10  A  cubic  box.  As  shown  in  Fig.  lc,  the  calculated 
migration  energy  of  LSM50  is  0.76  eV.  If  predicted  oxygen-vacancy 
formation  energy  of  3.25  eV  is  added,  it  is  closer  to  the  experimental 
result  (4.01  eV  in  Table  1).  Thus,  we  assumed  that  our  compu¬ 
tational  approach  using  the  surface  model  (see  Fig.  la)  may  be 
sufficient  enough  for  further  DFT  calculations.  After  replacing  Mn 
with  Cr,  Fe,  or  Co,  we  carried  out  similar  predictions  as  summarized 
in  Table  1 .  For  LSC50,  a  predicted  Ea  of  1 .87  eV  is  slightly  higher  than 
that  determined  from  experiment  ( 1 .41  ±  0.02  eV)  [12,13].  Our  pre¬ 
diction  clearly  shows  that  LSC50  has  the  lowest  diffusion  barrier  for 
oxygen  ion  transport  (LSC50  >  LSF50  >  LSM50  >  LSCr50),  suggesting 
that  LSC-based  MIEC  cathodes  allows  for  faster  ionic  transport  than 
other  three  materials. 

3.2.  Interactions  between  oxygen  species  and  Lao.5Sro.5BO2.75 

As  mentioned  earlier,  understanding  the  mechanisms  of  the 
interactions  between  oxygen  species  and  MIEC  cathode  surfaces 
are  crucial  to  achieving  rational  design  of  high-performance  cath¬ 
ode  materials  for  SOFCs.  As  Greeley  et  al.  [39]  described  on  the  basis 
of  the  Sabatier  principle  [40],  the  bond  strength  between  interme¬ 
diate  species  and  surfaces  plays  a  significant  role  in  understanding 
catalytic  activity  of  surfaces.  To  examine  the  interactions  of  oxy¬ 
gen  species  with  four  MIEC  cathode  materials,  we  constructed  the 
(10  0)  surface  for  each,  which  is  energetically  the  most  stable  one 
among  the  low-index  (1 1  0),  (1 1  0),  and  (1  0  0)  surfaces  [8,23]. 
Then,  we  carried  out  2D  slab-model  calculations  based  on  the  four- 
layer  (10  0)  surfaces  to  examine  O2-Lao.5Sro.5BO2.75  interactions, 
as  shown  in  Fig.  2a.  A  vacuum  space  with  ~24  A  in  the  perpendic¬ 
ular  direction  to  the  surface  was  used  to  guarantee  no  interactions 
between  slabs.  All  surface  calculations  were  performed  by  relaxing 
the  top  two  layers  while  fixing  the  bottom  two  layers  to  the  esti¬ 
mated  bulk  parameters.  The  relative  energies  for  the  02 -surface 
interactions  are  calculated  according  to  AEads  =  EP  -ER,  where  EP 
and  Er  are  the  calculated  energies  of  products  and  reactants,  respec¬ 
tively.  Illustrated  in  Fig.  2b  are  energetically  the  most  plausible 
configurations,  including  molecular  adsorption,  dissociation,  and 
incorporation  into  an  oxygen  vacancy.  Table  1  compiles  the  adsorp¬ 
tion  energies  of  energetically  stable  configurations  for  kinetic  and 
mechanistic  studies  [9].  Further,  Fig.  2c  clearly  shows  that  the 
adsorption  energies  of  adsorbed  species  (e.g.,  molecular  adsorption, 
dissociation,  and  incorporation  processes)  are  well  correlated  with 
the  B  cations  in  the  four  Sr-doped  LaB03  cathodes.  This  means  that 


the  bond  strength  between  B  cation  and  0  atom  may  be  a  crucial 
indication  of  catalytic  activity  toward  oxygen  reduction  for  charge 
transfer  from  d  orbital  of  B  cations  to  2p  orbital  of  molecular  oxygen 
[7].  Cr-based  MIEC  cathodes  have  the  strongest  bond  with  adsorbed 
oxygen  species.  However,  regarding  surface  exchange  coefficients 
associated  with  a  surface  catalytic  activity  toward  oxygen  reduc¬ 
tion,  LSC50  has  a  much  faster  02  kinetics  than  LSM50  [12,13]. 
Although  one  has  to  take  into  account  other  parameters  (e.g.,  the 
concentrations  of  oxygen  vacancies,  electrons,  and  electron  holes) 
[41  ]  to  estimate  surface  oxygen  exchange  coefficients,  the  trend 
shown  in  Fig.  2c  seems  to  suggest  that  the  adsorption-energy  cal¬ 
culations  are  reasonably  related  to  the  catalytic  activity  toward 
oxygen  reduction.  Compared  with  the  parameters  determined  from 
SIMS  [12,13],  our  DFT  calculations  verify  that  LSC-based  MIEC 
cathodes  have  the  highest  catalytic  activity.  Furthermore,  Fig.  2d 
manifests  that  the  surface  diffusion  of  adsorbed  oxygen  species 
after  dissociation  also  have  a  strong  relevance  to  the  B  cations,  indi¬ 
cating  that  the  bond  strength  between  surfaces  and  oxygen  species 
becomes  weaker  as  Z  increases.  This  may  mean  that  the  adsorbed 
oxygen  species  or  oxygen  ions  after  dissociation/incorporation  pro¬ 
cesses,  respectively,  on  LSC50  can  more  easily  diffuse  on  the  surface 
or  to  the  TPB  compared  to  others.  In  addition,  on  the  basis  of  the  esti¬ 
mation  of  the  bulk-diffusion  barriers  (Ea)  of  oxygen  ions  through 
the  MIEC  bulk  phases  and  the  adsorption  energies  (AEads  02), 
we  examined  their  correlation  as  shown  in  Fig.  3.  The  four 
cathodes  materials  (Lao.5Sro.5CrO2.75,  LSCr50;  Lao.5Sro.5MnO2.75, 
LSM50;  Lao.5Sro.5FeO2.75,  LSF50;  Lao.5Sro.5CoO2.75,  LSC50)  show  a 
very  good  correlation  with  a  slight  deviation  of  LSM50. 

We  then  examined  the  effect  of  A  cation  substitution  by 
replacing  La  with  Ce,  Pr,  Nd,  Pm,  and  Sm  (Ceo.5Sro.5MnO2.75, 
CeSM50;  Pro.5Sro.5MnO2.75,  PrSM50;  Ndo.5Sro.5MnO2.75,  NdSM50; 
Pmo.5Sro.5MnO2.75,  PmSM50;  Smo.5Sro.5MnO2.75,  SSM50,  respec¬ 
tively).  Interestingly,  only  SSM50  have  a  much  closer  character  to 
LSC50,  implying  that  Sm  may  play  a  role  in  exhibiting  a  higher 
catalytic  activity  toward  oxygen  reduction  and  a  lower  diffusion 
barrier.  Our  further  examination  of  Smo.5Sro.5CoO2.75  (SSC50)  and 
Smo.5Sro.5FeO2.75  (SSF50)  verifies  the  Sm  effect  (see  Fig.  3).  It  is 
noted  that  our  prediction  of  the  activation  barrier  of  SSC50  is  in 
excellent  agreement  with  the  experimental  value  (1.1  eVvs.  0.9  eV) 
[42].  Thus,  Fig.  3  manifests  that  the  adsorption  energies  of  oxy¬ 
gen  species  have  a  strong  correlation  with  the  diffusion  barriers  of 
oxygen  ion  transport.  Comparison  with  cathode  materials  that  are 


Fig.  3.  Comparison  of  adsorption  energies  versus  diffusion  barriers  of  oxygen  ions 
through  the  MIEC  bulk  phases.  Ea  denotes  the  diffusion  barriers  of  oxygen  ion 
transport  through  the  bulk  MIEC  phases,  while  AEads  02  is  the  averaged  adsorp¬ 
tion  energies  of  superoxo-  and  peroxo-like  species  on  B  cations.  The  linear  line  in 
red  was  obtained  from  LSCr50,  LSM50,  LSF50,  and  LSC50. 
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known  catalytically  active  at  low  temperatures  (LSC  and  SSC)  sug¬ 
gests  that  tailoring  novel  cathode  materials  must  also  consider  the 
need  for  a  weak  binding  between  a  surface  and  oxygen  species  with 
a  low  O-diffusion  barrier,  although  other  factors  (such  as  thermal 
mismatch  with  the  electrolyte  materials)  have  to  be  considered  as 
well.  We  will  extend  our  simulations  to  more  realistic  surface  mod¬ 
els  (Ai _xA^)(Bi_yBy)03_5  (i.e.,  A  =  La  or  Sm,  A'  =  Sr,  B  =  Co,  Br  =  Fe  or 
Mn)  to  predict  the  effect  of  defects  and  other  relevant  conditions 
on  surface  catalytic  and  ionic  transport  properties. 

4.  Conclusion 

Our  comparative  study  demonstrated  that  LSC  is  the  most  effi¬ 
cient  cathode  material  among  the  four  considered  due  to  its  fast 
diffusion  on  surfaces  and  through  the  bulk  phases,  qualitatively 
in  agreement  with  experimental  results  derived  from  the  perfor¬ 
mances  of  SOFCs  based  on  the  materials  [6].  Prediction  of  diffusion 
barriers  and  adsorption  energies  after  dissociation/incorporation 
may  provide  vital  information  for  rational  design  of  better  ABO3- 
type  cathode  materials  for  SOFCs.  Thus,  DFT  calculations  can  be  a 
useful  tool  to  search  for  high-performance  cathode  materials.  Fur¬ 
ther,  DFT  based  results  may  be  combined  with  kinetic  Monte  Carlo 
(KMC)  methods  to  simulate  electrochemical  behavior  of  cathode 
materials  under  SOFC  operation  conditions  [33,43],  which  can  be 
verified  by  electrochemical  measurements. 
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